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Abstract. Differential properties for orthogonal polynomials in several vari- 
ables are studied. We consider multivariate orthogonal polynomials whose 
gradients satisfy some quasi-orthogonality conditions. We obtain several char- 
acterizations for these polynomials including the analogous of the semiclas- 
sical Pearson differential equation, the structure relation and a differential- 
difference equation. 



1. Introduction 

In the univariate semiclassical moment functional u is a quasi-definite 

functional satisfying a distributional Pearson equation 

(1) D{(j)u) = ipu 

where <fr and ip are polynomials with deg(^) > 1. They constitute a natural ex- 
tension of classical functionals (Hermite, Laguerre, Jacobi, and Bessel) and they 
have been extensively analyzed during the last two decades (see [H], Obvi- 
ously, the sequence {p n } n of orthogonal polynomials associated with u is also called 
semiclassical. 

The distributional Pearson equation 0) plays a key role in the study of dif- 
ferential properties of semiclassical orthogonal polynomials. In fact, this equation 
allows us to characterize semiclassical polynomials as the only sequences of orthog- 
onal polynomials satisfying one of the following equivalent properties: 

(a) the so-called structure relation, 

(b) the quasi-orthogonality of the derivatives, 

(c) the second order differential-difference relation. 

For further properties and characterizations see, for instance, |15j . 

In the multivariate case, we will call semiclassical a quasi-definite moment func- 
tional u satisfying the matrix Pearson-type equation 

div ($ u) = \ff* u, 

where $ is a d x al symmetric polynomial matrix and W is a d x 1 polynomial matrix 
with deg^ > 1, and such that det (it, $} ^ 0. Of course, this definition includes all 
the "classical" moment functionals in the usual literature ( |21 151151 H0H13lH^lll6| ). 
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The aim of our contribution is to analyze the extension of characterizations (a), 
(b), and (c) to the multivariate case by means of a matrix formalism. In the 
bivariate case, the hrst of these three characterizations was the main result of a 
previous publication (see £Q). The two other characterizations constitute the main 
objective of this work. 

The structure of the paper is the following. In Section 2, a basic background 
about moment functionals and multivariate orthogonal polynomials is given in or- 
der to allow the reader to be familiar with such concepts. In Section 3 we present 
two new characterizations for sequences of multivariate semiclassical orthogonal 
polynomials. Theorems 13.31 and 13.41 are, respectively, the extensions of the quasi- 
orthogonality of the derivatives and the differential-difference relation satisfied by 
semiclassical orthogonal polynomials in one variable. Finally, in Section 4 we ex- 
plore two non-trivial examples of multivariate semiclassical orthogonal polynomials. 



2. Background 

We denote by V the linear space of polynomials in d variables with real coeffi- 
cients, and by V its topological dual (see ^H]). Let V n be the set of real polynomials 
of total degree not greater than n, and H n the space of homogeneous polynomials 
of degree n in d variables. It can be observed that 

dimT 5 ,! =(^| and dimLI n = ( 

\ n ) \ n 

Let Aihxki^-) be the linear space of (h x k) real matrices, and A / J/ lX fc('P) the linear 
space of (h x k) polynomial matrices. In addition, 1^ represent the identity matrix 
of order h. 

If M — {fni,j(x,y)) i ' : j_ 1 represents a (h x k) polynomial matrix, we define the 
degree of M by 

degM = max{degmij(x,y), 1 < i < h, 1 < j < k} > 0. 

Now, we review some basic definitions and tools about multivariate orthogonal 
polynomials that we will need in the rest of the paper. For a complete description 
of this and other related subjects see, for instance, [Tl Efl 151 ITT1 1121 1161 117j . 

Let No be the set of nonnegative integers. For a = (oti, . . . ,atd) G Ng and 
x = (xi, . . . , Xd) € M d we write x" = a;" 1 ■ • ■ x^ d . The number |a| = ai + • • • + ay 
is called the total degree of x a . 

Let {^a} Q £Ng be a multi-sequence of real numbers and let us denote by u, the 
only linear functional in V' satisfying 

(u,x a ) = fl a . 

where, as usual, ( , } stands for the duality bracket. Then, u is called the moment 
functional determined by {^ a } a &i^ an< ^ ^ ne number [i a is called the moment of 
order a. 

Let (P™)| a | =n be a sequence of polynomials of total degree n. Using the matrix 
notation introduced in |lllll2j and [Uj, we can denote by V n the vector polynomial 

P„ = ( J P Q ")| Q |=n = (O) . O) ' * ■ ■ ' P a(^Y € M rn xl(Pn), 

where a' 1 ', a^ 2 \ ■ • • , a' r "^ are the elements in {a £ Nq : \a\ = n} arranged accord- 
ing to the lexicographical order. 
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When {P m }m=o i s a basis of V n for each n > 0, then {P n }„>o is called a 
polynomial system (PS). 

Definition 2.1. We will say that a PS {P n } n >o is a weak orthogonal polynomial 
system (WOPS) with respect to u if 

(] ( M ,P„P*„) = 0, m^n, 

where H n 6 Al r „(M) is a non singular matrix. 

This definition means that every polynomial component in P„, is orthogonal to 
all the polynomials of lower degree, but two polynomial components of the same 
degree don't have to be orthogonal. 

A moment functional u is said quasi-definite if there exists a WOPS with respect 
to u (0). 

Now, we are going to recover the three term recurrence relation for orthogonal 
polynomials in several variables (pi lllllT^lll7| ). This relation takes a vector- matrix 
form and it plays an essential role in understanding the structure of orthogonal 
polynomials, as in the univariate case. 

For n > 0, there exist unique matrices A n ^ £ A4 r „ X r„ + i(K), B n ,i G _M r „ X r„(R) 
and C nt i G Mr n xr n ^ 1 (K), i — 1> 2, . . . , d, such that 



(3) Xi P n — A n> i P n +i + B n ^ P n + C n> i P„-i, 



1 < i < d. 



where P_i = and C_i ;i = 
rank A 



Moreover, rank^4„i = rankC* 



where A r . 



rankC* +1 = r n+1 , 



and C£ +1 



V ^n+l,d 

As a consequence of the three term recurrence relation, it is possible to obtain 
the polynomial P„+i in terms of the vector polynomials P„ and P n _i (see [H]). 

Let = (D^ x . . . £)* d ) G A4 rn+1 xdr„(R) be a left generalized inverse of A n . 
There exist matrices E" +1 and E%^\ such that 



(4) 



d 



- 71+1 



/ / •£» D n j P n 



n+1 ■ 



■ E n+1 



_1 Jr«-1- 



Moreover, since the matrices C n j, for i = 1, 2, . . . , d, are full rank, we can obtain 



P ra _i in terms of the vector polynomials P„+i and 
generalized inverse of C n ,i, that we will denote by G 
relations: 



n • 



In fact, using the left 
we deduce the following 



(5) P„-i = f P n+1 + ( Xi G n ,i - G nti B n ,i) P„, 1 < i < d. 

3. SEMICLASSICAL MULTIVARIATE ORTHOGONAL POLYNOMIALS 

First, we introduce the concept of semiclassical moment functional. For d = 2, 
this definition was given in pp. 
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Definition 3.1. A quasi-definite moment functional u is said to be semiclassical 
if it satisfies the matrix Pearson-type equation 

(6) div ($ u) = u, 
where 

* = (<fe)- J= i G Md ^> * = (^)ti € M dx i(P), 
are polynomial matrices such that $ is symmetric, deg<J> = p > 0, deg 'I' = g > 1, 
and 

(7) det(u,$}^0. 
We denote by s = max{p — 2, q — 1} > 0. 

Expression means 

(div ($ «),/) = <** u,/>, 

that is, 

<u, $ v/ + * /) = o, v/ e p. 

The natural extension of the above property for matrices involves the Kronecker 
product (see, for instance [7], p. 242), 



(8) div (($ ® J h ) u) = ® 4) u, fc>l. 
Relation |JS} is equivalent to 

(9) ($ (8 7 h ) V u = ® I h ) u, h>l, 
where $ = $ — (div $)*. 



Remark If s = 0, that is, deg $ = p < 2 and deg $ = 1, we recover the definition 
of classical WOPS given in 0, which includes the classical bivariate orthogonal 
polynomials studied by H. L. Krall and I. M. Sheffer (^31), and other authors 

(dm mi mum). 

A WOPS with respect to a semiclassical moment functional u is called semiclas- 
sical. 

Now, we are going to prove three characterizations for multivariate semiclassical 
orthogonal polynomials: structure relation, quasi-orthogonality relation for the 
gradients, and a differential-difference relation. From now on, we will denote by 
{I\i}n>o a given WOPS associated with a quasi-definite moment functional u. 

Theorem 3.1 (Structure relation). The moment functional u is semiclassical if 
and only «/{P,i}ri>o satisfy 

n+p— 1 

(10) $ VPj, = {Id® PpF™, for n > s + 1, 

j=n— a-l 

where Ff e M drj xrS^)- 

Proof. See Theorem 5 in pQ. □ 

Using relations (@J and 10, a shorter structure relation whose coefficients are 
polynomial matrices can be deduced. 

Corollary 3.2. If u is semiclassical, then {P n }n>o satisfy 

$ VP^ = (I d <g> Fi +1 )M? + (I d ® K)M£, for n > a + 1, 
where M™ are polynomial matrices with deg(M") < s and deg(Af^) < s + 1. 
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From Theorem 13. II we can obtain the second characterization for semiclassical 
orthogonal polynomials. 

Theorem 3.3 (Quasi-orthogonality relation for gradients). The functional u is 
semiclassical if and only if, for n > s + 1, 

(u, (VP*J* $ VP*,) = 0, < m < n - s. 

Proof. Suppose that u is semiclassical. From Hl()|l . we can write 

n+p—l 

(«,(VPt,)**VF*)= £ MVP*J*(I<i®P*)>i^. 

j—n—s—1 

Taking into account that 

(vpL) 4 (/d®P5) = ( (9iP m )P*- 1 (5 2 p m )p* I ••■ I (a d p m )p* ) , 

we obtain 

( lt ,(VP^) i (/ rf ®P*)) = 0, m-Kj, 

that is, 

(u, (VP*J*$VP*J = 0, m<n-s. 
Reciprocally, assume that the quasi-orthogonality relations hold, we define 

s+l 

* = -^2(u, (VP*)* $ VP*)^ 1 p s: . 

i=0 

Since VP^ — I2, we obtain 

(div ($ u),P*J = -(it,$VP£) = -(u,(Wl)'$VP£). 
On the other hand, 

s+l 

(** u,P*J = <u,*P^) = -^(^(VPi^VP^tfr^p^). 

i=0 

If n > s + 2, using the above relations, we get: 

(div (<& u),P£) = () = <** u,P*J. 
Furthermore, for < n < s+ 1, 

(div ($ «),!*) = -( U) (VPi)*$VP*,) = (** u,P^), 
and then, the result follows. □ 

Now, we deduce a matrix differential-difference relation for semiclassical orthog- 
onal polynomials in several variables. 
Let us define the differential operator 

L[f] = div ($V/) + tf*v/, VfeV. 

Therefore, the Lagrange adjoint of L is given by 

(11) L*[u] = div ($Vti) - div (*u), 

since it satisfies 

(L*[u]J) = (u,L[f}), VfGP. 
Using the explicit expression of the polynomial matrices $ and 'F, the operator 
L can be written as follows 
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d d 

L[f}= E • y, 1 ■■<>■■>'■ 

i,j=l i=l 

and, from the above expression, we deduce degL[/] < s + deg/. 

Theorem 3.4 (Matrix differential -difference relation). A functional u is semiclas- 
sical if and only if there exist matrices A™ £ A4 ri xr„(^), such that 

n+s 

(12) LK] = J2 P * A ?> n>s + l. 

i—n — s 

When n < s, relation %Tgj} reads 

n-\-s 
i=l 

that is, Ag = 0,Vn > 0. 

Proof. If u is a semiclassical functional, Lemma 4.f in 5 provides 

<u,P m div ($VP* t )) = <u,div ((I d ®P m ) <& VP* n )> - («, (Wj* $ VP* n >, 
for u, m > 0. Besides, from (|UJ|, we deduce 

(«, div ((7 d ® P m )$VP*J) = -<V«, ($ ® / rm )(/ d ® P m )VP*„) 

= -<($®I rm )Vu,(/ rf ®P m )VP^> 
= -<(*®/ rm )«,(/ (J ®P TO )VP*> 

= -<U,(**(»/r m )(id®Pm)VP^> 

= -<u,P m **VP^>. 

Therefore, we have 

(13) <^P m «) = -(«, (VP'„)'$ VP,\). 

Now, since £[P£J is a polynomial matrix of degree at most n + s, we obtain the 
expansion 

n+s 

i[p* l ] = E^Ar, 

i=0 

where AJ^ is given by 

n+s 

(u, P m L[P*J) = (u, P m £ F * A D = ^ m Ar„, 

z=0 

and using Q13|). we get 

F ro A^ = -< U ,(VP^) t $VP^). 
In the case n > s + 1, from Theorem l3.3l 

iJ m A™ = 0, < m < n - s, 

so, relation H12fl follows. 
When n < s, 

H Q AS = -(«,(VPS) t $ VP* t ) =0. 
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Reciprocally, if relations Ijl2|) and (|13|l hold, then, L*[u] = 0, since 

n+s 

(L*M,P* t ) = («,LpP*]) = < U > P *> A " = °' 

i-max{Ln-s} 

Therefore, Lemma 3.4 in gives 

= (L*[u],F m V t n ) = (u,L[V m F t n \) 

= (u, L[F m ] P<J + (u,P m L[1P* ]> + 2<u, (VPjJ* $ VP*„). 

If < m < n - s, (u,L[F m ] P*,} = 0, since deg(L[P m ]) <m + s. On the other 
hand, from (|12|) . we deduce (w,P m = 0, and therefore, we conclude 

(u,(W¥ t J t ^V¥ t n } = 0. 

Finally, Theorem 13 . 31 provides the desired result. □ 

Again, three term recurrence relations (@J and 0, allow us to express £[PJJ in 
terms of the vector polynomials P^ +1 and P^: 

Corollary 3.5. If a quasi-definite moment functional u is semiclassical, then there 
exist polynomial matrices N™ , and N£ , satisfying 

(14) L[P*]=P^ +1 iVf + P^iV 2 " 

with deg{N[ l ) < s - 1 and deg(A 2 ") < s. 

4. Examples 

Example 1: Appell type orthogonal polynomials 

The so-called Appell polynomials ([5]) are orthogonal polynomials in d variables 
associated with the weight function 

Wq (x)=x q (1-|x|)^ 

where a = (a%, . . . ,a d ) G R d , with cti > —1, 1 < i < d, and /3 G IR such that 
[3 > — 1, on the cZ-simplex, 

T d = {x = (xi, . . . ,x d ) : xi>0,...,Xd>0, 1 - |x| > 0}. 

The Appell moment functional u is defined as follows 



(«,/> 



/(x) w Q (x)dx. 



Using our definitions (PQ,[S]), u is a classical moment functional (i.e., it is semiclas- 
sical with s = 0) , since it satisfies the matrix Pearson-type equation JSJ , where the 
matrices $ and "J are given by 

/xi(a;i-l) xix 2 ••• xix d \ 
x 2 x 1 x 2 (x 2 -l) ■■■ x 2 x d 



X d X 2 



X d {x d - 1)/ 



* = 



f{\a\ + d)x 1 - (an + \)\ 
'\a\ + d)x 2 - (a 2 + 1 



\(\a\+d)x d - {a d + l)J 
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Now, we introduce the Appell-type polynomials as the orthogonal polynomials in d 
variables with respect to the moment functional 

v = u + A <5(x), 

where A > is a positive real number, and S(x) is the usual Dirac distribution at 
OeK d . The action of v over polynomials is defined as follows, 

(«,/>= f /(x)w„(x)d* + A/(0). 

The moment functional v is semiclassical with s = 1, since v satisfies the matrix 
Pearson-type equation @ 

(15) div ($ v) = ¥ v, 

where 

$ = xi$, * = - l),xix 2 ,- • • ,xix d )* + xi*. 

In fact, using that xi 5(x) = 0, we get 

div (<l u) = div (xi $ (u + A5(x)) = 

= div (xi & u) + div (xi $ A <5(x)) = 
= (1,0, ••• ,0)$ M + xidiv ($ u) = 

= (xi(xi - 1), XiX 2 , • • • ,XiXd)v + Xi u. 

Observe that the matrix Pearson-type equation for i> is not unique. The moment 
functional v also satisfies (|15[) with 

$ = X,$, # = (xiXi, • • • ,Xi(Xi - 1), • • • ,XjXd)* + Xj#, 

for 1 < i < d, since Xj <5(x) = 0. 

Example 2: A multivariate analogue of the classical orthogonal polyno- 
mials 

Examples of two-variables analogues of the Jacobi polynomials are studied in ^H] 
by T. Koornwinder. Using similar tools, we present an example of a semiclassical 
weight function with unbounded support. 

Let a, be real numbers with at > —1, for 1 < i < d. Then, for fcj > 0, 1 < i < d, 
and fci > kd, we define the polynomials 

p {a u ...,a d ), \ T { ai +2k d + l), ir(aj), s r (a«,_ 1 ), n fc d p (a dl 0) f_-l„,_N 

where £jj. a *^ is a Laguerre polynomial in one variable, and P^ d ^ is a Jacobi poly- 
nomial. 

The polynomials P^'k"^^) are orthogonal with respect to the weight function 
w(x) = x? 1 . . . x a d d _-i e^ {xi+ - +Xd - l) (1 - x^x d ) aa , 

on the region {x = (xi, . . . , xa) / — x\ < Xd < x\, Xj > 0, i = 1, . . . , d — 1}. 
Defining the matrices 

/ xi(xi - xa) O \ 

x 2 

Xd-l 

\ O x i(^i - x d ) ) 
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(-x\ + xi x d + (at + 2)xi + (a d 



ai 



l)x d \ 



Ot2 - X 2 



ad-i - 
-(a<j 



Xd-l 

-l)x\ 



we can prove that the weight function w(x) satisfies © and so, the polynomials 
fjt™ 1 '"^' (x) are semiclassical. 

Another examples for semiclassical orthogonal polynomials in two variables ap- 
pear in pp. 
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